function [varargout]=Equ_NK(PP,SS,order)

%'Eta_Z','Eta_Zp','Eta_M','Eta_Mp','N','w','mc','Y',...
 %'Yp','Profit','PAC','L','C','Pir','Pirp','T',...
 %'TAU','B','Lag_B','Lag_Bp','ir','Lag_ir','Lag_irp','rr',...
 %
 
 PP_RHO_Z=PP.RHO_Z;
PP_RHO_M=PP.RHO_M;
PP_ALPHA=PP.ALPHA;
PP_EPSILON=PP.EPSILON;
PP_THETA=PP.THETA;
PP_Taylor_Pi=PP.Taylor_Pi;
PP_Taylor_ir=PP.Taylor_ir;

 SS_Z=SS.Z;
SS_N=SS.N;
SS_w=SS.w;
SS_mc=SS.mc;
SS_Y=SS.Y;
SS_Profit=SS.Profit;
SS_PAC=SS.PAC;
SS_C=SS.C;
SS_L=SS.L;
SS_Pir=SS.Pir;
SS_ir=SS.ir;
SS_TAU=SS.TAU;
SS_B=SS.B;
SS_T=SS.T;
SS_G=SS.G;

 Eta_Z=0;
Eta_Zp=0;
Eta_M=0;
Eta_Mp=0;
N=0;
w=0;
mc=0;
Y=0;
Yp=0;
Profit=0;
PAC=0;
L=0;
C=0;
Pir=0;
Pirp=0;
T=0;
TAU=0;
B=0;
Lag_B=0;
Lag_Bp=0;
ir=0;
Lag_ir=0;
Lag_irp=0;
rr=0;
if order==0
res=ones(15,1);
res(1)=Eta_Z*PP_RHO_Z - Eta_Zp;
res(2)=Eta_M*PP_RHO_M - Eta_Mp;
res(3)=SS_Z*exp(Eta_Z)*(SS_N*exp(N))^(1 - PP_ALPHA) - SS_Y*exp(Y);
res(4)=- SS_mc*exp(mc) - (SS_N*SS_w*exp(-Y)*exp(N + w))/(SS_Y*(PP_ALPHA - 1));
res(5)=(Pirp*(SS_Pir + 1))/(SS_ir + 1) - Pir + (PP_EPSILON*mc)/(PP_THETA*(SS_Pir + 1));
res(6)=-PAC;
res(7)=SS_Y*exp(Y) - Profit - SS_PAC - SS_Profit - PAC - SS_N*SS_w*exp(N + w);
res(8)=Eta_M - ir + PP_Taylor_Pi*Pir + Lag_ir*PP_Taylor_ir;
res(9)=(B + SS_B)/(SS_ir + ir + 1) - SS_T - T - (Lag_B + SS_B)/(Pir + SS_Pir + 1) - SS_G + SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);
res(10)=-B;
res(11)=-T;
res(12)=ir - Pirp - rr;
res(13)=SS_L*exp(L) - SS_N*exp(N);
res(14)=B - Lag_Bp;
res(15)=ir - Lag_irp;
varargout{1}=res;
elseif order==1
fir_ord.Eta_Z=[PP_RHO_Z;0;SS_Z*exp(Eta_Z)*(SS_N*exp(N))^(1 - PP_ALPHA);0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Zp=[-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M=[0;PP_RHO_M;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.Eta_Mp=[0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N=[0;0;-(SS_N*SS_Z*exp(Eta_Z)*exp(N)*(PP_ALPHA - 1))/(SS_N*exp(N))^PP_ALPHA;-(SS_N*SS_w*exp(-Y)*exp(N + w))/(SS_Y*(PP_ALPHA - 1));0;0;-SS_N*SS_w*exp(N + w);0;0;0;0;0;-SS_N*exp(N);0;0];
fir_ord.w=[0;0;0;-(SS_N*SS_w*exp(-Y)*exp(N + w))/(SS_Y*(PP_ALPHA - 1));0;0;-SS_N*SS_w*exp(N + w);0;SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);0;0;0;0;0;0];
fir_ord.mc=[0;0;0;-SS_mc*exp(mc);PP_EPSILON/(PP_THETA*(SS_Pir + 1));0;0;0;0;0;0;0;0;0;0];
fir_ord.Y=[0;0;-SS_Y*exp(Y);(SS_N*SS_w*exp(-Y)*exp(N + w))/(SS_Y*(PP_ALPHA - 1));0;0;SS_Y*exp(Y);0;0;0;0;0;0;0;0];
fir_ord.Yp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Profit=[0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0];
fir_ord.PAC=[0;0;0;0;0;-1;-1;0;0;0;0;0;0;0;0];
fir_ord.L=[0;0;0;0;0;0;0;0;SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);0;0;0;SS_L*exp(L);0;0];
fir_ord.C=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir=[0;0;0;0;-1;0;0;PP_Taylor_Pi;(Lag_B + SS_B)/(Pir + SS_Pir + 1)^2;0;0;0;0;0;0];
fir_ord.Pirp=[0;0;0;0;(SS_Pir + 1)/(SS_ir + 1);0;0;0;0;0;0;-1;0;0;0];
fir_ord.T=[0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0];
fir_ord.TAU=[0;0;0;0;0;0;0;0;SS_L*SS_w*exp(L + w);0;0;0;0;0;0];
fir_ord.B=[0;0;0;0;0;0;0;0;1/(SS_ir + ir + 1);-1;0;0;0;1;0];
fir_ord.Lag_B=[0;0;0;0;0;0;0;0;-1/(Pir + SS_Pir + 1);0;0;0;0;0;0];
fir_ord.Lag_Bp=[0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0];
fir_ord.ir=[0;0;0;0;0;0;0;-1;-(B + SS_B)/(SS_ir + ir + 1)^2;0;0;1;0;0;1];
fir_ord.Lag_ir=[0;0;0;0;0;0;0;PP_Taylor_ir;0;0;0;0;0;0;0];
fir_ord.Lag_irp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1];
fir_ord.rr=[0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0];
varargout{1}=fir_ord;

 end
